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We outline the proofs of the factorization theorems for exclusive two-body charmonium production 
in B-meson decay and e^e~ annihilation to all orders in perturbation theory in quantum chromo- 
dynamics. We find that factorized expressions hold up to corrections of order mc/rub in _B-meson 
decay and corrections of order rn^/s in e^e~ annihilation, where nic is the charm-quark mass, mi, 
is the bottom-quark mass, and ^/s is the e'^e~ center-of-momentum energy. 
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The exclusive decays of B mesons into a light meson 
plus a charmonium are significant tools for understanding 
the weak interactions. For example, they could provide 
new constraints on the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix and enhance our understanding of the 
origins of CP violation. However, the effects of the 
strong interactions are significant in such processes and 
present an obstacle to achieving reliable theoretical cal- 
culations of the process rates. The exclusive production 
of double-charmonium in e^e~ annihilation provides an 
arena in which to explore the mechanisms of charmo- 
nium production and the interplay between the pertur- 
bative and nonperturbative regimes of quantum chromo- 
dynamics (QCD). However, measurements of the cross 
sections for double-charmonium production by the Belle 
[H and BABAR [2] Collaborations have spurred a re- 
examination of the theoretical predictions for those cross 
sections 0. 

In theoretical computations of the rates for the exclu- 
sive decays of B mesons into a light meson plus a char- 
monium and of the rates for the exclusive production of 
double-charmonium in e"*'e~ annihilation, a crucial step 
is the separation of the effects of the strong interactions 
into short-distance, perturbatively calculable contribu- 
tions and long-distance, inherently nonperturbative con- 
tributions. Such separations are usually embodied in fac- 
torization theorems. In the case of the exclusive decays of 
B mesons into a light meson plus a charmonium, several 
factorization theorems have been posited 0, d, [1] . In the 
case of the exclusive production of double-charmonium 
in e^e~ annihilation, factorization conjectures have gen- 
erally been formulated in terms of nonrelativistic QCD 
(NRQCD) (3. In this letter, we outline the proofs of the 
factorization theorems for these processes in QCD to all 
orders in perturbation theory. 

We wish to show that the amplitudes for these pro- 
cesses can be written in a factorized form. The mean- 
ing of this statement, for the case of exclusive decays 
of B mesons into a light meson plus a charmonium, is 
that the decay amplitude is decomposed into the sum 
of products of a i?-meson-to-light-meson form factor and 



an amplitude for a charm-quark-antiquark (cc) pair to 
be produced at short-distances in a color-singlet state 
and evolve into a charmonium. The form factor con- 
tains a term that can be decomposed into a convolution 
of a hard-scattering amplitude with i3-meson and light- 
meson light-cone amplitudes and a term that cannot be 
decomposed further. (These two terms are analogous to 
the two terms in the factorization formula in Eq. (4) 
of Ref. Q for the case of decays to two light mesons.) 
The cc-to-charmonium amplitudes can be further decom- 
posed into a sum of products of long-distance NRQCD 
matrix elements times short-distance coefficients. We ar- 
gue that this factorized form holds up to corrections of 
relative order nic/y/s, where s = m^. This result was 
suggested previously in Ref. 01 ■ For the case of the 
exclusive production of double-charmonium in e~^e~ an- 
nihilation, factorization means that the production am- 
plitude decomposes into a sum of products of a short- 
distance amplitude and two amplitudes for a cc pair to 
be produced at short-distances in a color-singlet state and 
evolve into a charmonium. Again, the cc-to-charmonium 
amplitudes can be further decomposed into a sum of 
products of long-distance NRQCD matrix elements times 
short-distance coefficients. We argue that this factorized 
form holds up to corrections of relative order m'^/s, where 
is the center-of-momentum (CM) energy of the e"'"e~ 
pair. Although our analyses are for the specific cases of 
B decay and e"*"e~ annihilation, the techniques that we 
describe should apply to other exclusive quarkonium pro- 
duction processes and may also shed light on factoriza- 
tion in inclusive quarkonium production. We note that, 
because we consider exclusive two-body decays, we avoid 
the issues raised in Ref. 0] concerning light particles co- 
moving with the charmonium and the issues raised in 
Ref. concerning an additional heavy quark co-moving 
with the charmonium. 

We carry out our analyses in the _B-meson rest frame 
and in the CM frame of the e+e^ pair, choosing the 
three-momentum of the light meson or one of the char- 
monia to be in the negative 3 direction and choosing 
the three-momentum of the other charmonium to be in 
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the positive 3 direction. For a momentum k, we define 
light-cone momentum components = ± k^)/^/2, 
k±. We model the B meson as an on-shell active bot- 
tom quark, which participates in the electroweak inter- 
action, and an on-shell spectator light antiquark, which 
does not participate in the electroweak interaction. We 
take the quark and antiquark to be in a color-singlet 
state. We take the bottom quark to have momentum pf,, 
with = mb and all other components of momentum 
zero. We take the spectator quark to have momentum 
pi, with pz^-^Aqcd, the QCD scale. We neglect Aqcd 
and the difference between uib and the i?-meson mass 
in comparison with nib- Similarly, we model the light 
meson as an on-shell active light quark and an on-shell 
spectator light antiquark, with the quark and antiquark 
in a color-singlet state. We take the active-quark mo- 
mentum to be = ypk + qk and the spectator-quark 
momentum to be pk2 — (l~y)Pfc ~ <Zfc, where pk is the 
light-meson momentum, ^ mi,, = 0, pki_ = 0, 
qk± Aqcd, and pi ^ ^ K^^^-^/rnb. Finally, we 
model a charmonium as an on-shell cc pair in a color- 
singlet state, with the momentum of the c (c) equal to 
Pi {Pi)- We take pi = Pi/2 -\- qi and pi = Pi/2 - q^, 
where Pi is the charmonium momentum and Pi ■ qi — 0. 
In the charmonium rest frame, qi has only spatial com- 
ponents, whose magnitudes are of order rricV, where v 
is the typical charm-quark velocity in the charmonium 
rest frame (w^ « 0.3). In the e+e^ CM frame or B- 
meson rest frame, P^ ^ P2 ~ Q, Pii_ — P2± — 0, 



and Pf - P+ . 
momentum scale 



m1/Q, where Q = -^i is the large 
One can think of the on-shell ampli- 
tude that we use in our model as the on-shell perturba- 
tive QCD amplitude that is matched to a soft-coUinear- 
effective-theory amplitude in the case of the light meson, 
a heavy-quark-effective-theory amplitude in the case of 
the B meson, and an NRQCD amplitude in the case of 
the charmonia. 

The next step is to identify the regions of loop mo- 
menta that give contributions that are leading in powers 
of Q when we dress the lowest-order decay and produc- 
tion amplitudes in these models with additional gluons. 
These leading momentum regions correspond to partic- 
ular Feynman-diagram topologies [lo|, which we de- 
scribe below. The identification of leading regions has 
been discussed in detail in Ref. For a loop momen- 



tum k, the leading regions are the hard region, in which 
fc^ ~ Q; the collinear-to-plus (minus) region, in which 
|fc+|>|fc-| (|fc~|>|fc+|), k+k-r^kl<t:Q'^, and the soft 
region, in which k^ ^ k^ ~|fc±| <C Q- The collinear re- 
gions correspond to the directions of the momenta of the 
final-state light meson and charmonia. We note that, in 
the case of i3-meson decays, there is also a "semi-hard 
region" in which propagators are off shell by an amount 
of order mbAqcD- We treat this semi- hard region as part 
of the hard region. The "Glauber" region is also leading 
in power counting 12, In this region, |fc+| <ti \k±\. 





FIG. 1: Leading regions for the decay of a B meson to a light 
meson and a charmonium. 




FIG. 2: Leading regions for double-charmonium production 
in e'^e~ annihilation. The wavy line represents the virtual 



photon. 



^ fc^ ^ Q'^. However, for the exclusive pro- 
cesses we are considering here, the contours of integration 
of loop momenta are not pinched in the Glauber region, 
and it is possible to deform them out of it on a diagram- 
by-diagram basis. This is in contrast with the situation 
in, for instance, the Droll- Yan process, for which such a 
dia gram - by-diagram contour deformation is not possible 
[nHH,!!!- Therefore, we ignore the Glauber region in 
the remainder of our discussion. Contributions from the 
hard (semi-hard) region involve propagator denominators 
that are of order s (mbAqcD) and can be calculated in 
perturbation theory. Integrations of gluon momenta over 
the soft region or a collinear region that is associated with 
a massless particle encounter singularities in propaga- 
tor denominators that result in logarithmic divergences. 
[The soft singularities are at fc^ = 0, and the collinear 
singularities are at fc^ = 0, fej^ = (k^ = Q, k± = 0) for 
the collinear-to-plus (to-minus) region.] Thus, contribu- 
tions from the soft region and the collinear regions that 
are associated with a massless particle cannot be com- 
puted reliably in perturbation theory. The essence of 
the proofs of the factorization theorems that we consider 
here is to show that these contributions cancel or can be 
absorbed into the i?-meson-to-light-meson form factor or 
the nonperturbative NRQCD matrix elements. When 
a gluon three-momentum becomes parallel to the three- 
momentum of a charm quark or charm antiquark, the 
potential collinear singularity is shielded by rUc- Hence, 
the contributions from such collinear regions can be com- 
puted reliably in perturbation theory. Nevertheless, it is 
important for the arguments that we will make below to 
treat these regions separately from H. As we shall see, 
the contributions from these regions cancel. 

Now let us specify the diagrammatic topologies that 
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correspond to the leading regions. We work in the Feyn- 
man gauge. There are two distinct topologies in the 
case of B-meson decays: one in which the B-meson and 
light-meson spectators participate in the hard interac- 
tion and another in which they do not. These topologies 
are shown in Figs.[T]Ja) and[T{b), respectively. After one 
has achieved factorization, the contribution from the first 
topology can be decomposed into an expression that con- 
tains the convolution of a hard-scattering amplitude with 
_B-meson and light-meson light-cone amplitudes, while 
the contribution from the second topology cannot be de- 
composed further. The topology for the case of double- 
charmonium production is represented in Fig. [21 The 
components of the topologies are as follows: a B me- 
son (B); jet sub-diagrams for each of the coUinear re- 
gions, corresponding to a charmonium (Ji) and a light 
meson ( J2) in Fig. [T] and the two charmonia (Ji and 
J2) in Fig. [21 a hard sub-diagram {H) that includes the 
lowest-order annihilation or production process; and a 
soft sub-diagram (S). In H, all propagator denominators 
are of order s or TOhAgcD- Ji contains the active- and 
spectator-quark lines for a given meson or charmonium, 
as well as gluons and loops involving quarks and ghosts 
with momenta collinear to the meson or charmonium. Ji 
attaches to H through the active- and spectator-quark 
lines in the topology of Fig. [Ija), through the active- 
quark lines in the topology of Fig. [T{b) and through any 
number of gluons. S includes gluons with soft momenta 
and loops involving quarks and ghost with soft momenta. 
S attaches to Ji and to the B-meson quark and antiquark 
lines through any number of soft gluon lines. Recall 
that the _B-meson spectator antiquark itself carries a soft 
momentum. In Fig. [TJb), the light-antiquark-spectator 
line carries a momentum with minus component of order 
■Aqcd- Such a momentum can arise either through an 
endpoint contribution, in which p^T^ ~ Aqcd 0, [isl , or 
through the soft-coUinear-messenger mechanism |16| . in 
which an emitted gluon carries away most of ■ 

At this point we can outline the proof of factoriza- 
tion. Suppose that a gluon from Ji enters H. Then 
we can apply a collinear approximation to that gluon 
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I3L Il4l |. Specifically, we replace g^^^, in the gluon- 
propagator numerator with k^fii^/k ■ fii, where k is the 
gluon momentum, fii (712) is a unit light-like vector in 
the minus (plus) direction. The index fi corresponds to 
the attachment of the gluon to H , and the index v cor- 
responds to the attachment of the gluon to Ji. For k in 
the collinear region, the collinear approximation is valid 
at leading order in Q. Having made the collinear ap- 
proximation, one can use the factor k^ and the diagram- 
matic Ward identity (Feynman identity) for fermion lines 
^ ' 7 — -I- fc) • 7 -|- m] — [g • 7 -t- m] and its generalizations 
to non-Abelian gauge theories to show t hat the gluons 

nil 



with collinear momenta decouple from H [11 



In 

general, the decoupled gluons attach to light-like eikonal 
lines (path integrals of the gauge field) that connect to 



the point at which an active parton enters H. However, 
in our case these eikonal lines cancel because the light 
meson and the charmonia are color-singlet states. Now 
suppose that a gluon from S attaches to a jet sub-diagram 
at a line with momentum p. We can apply a soft approx- 
imation to that gluon 17[ • Specifically, we replace gp^i, in 
the gluon-propagator numerator with k^p^/k ■ p, where 
k is the gluon momentum, the index /x corresponds to 
the attachment of the gluon to Ji and the index v corre- 
sponds to the attachment of the gluon to S. We would 
like to use the factor fc^ in the soft approximation and the 
diagrammatic Ward identities to decouple the soft glu- 
ons. However, in general, the momentum p is different 
for every line to which the soft gluon attaches. That is, 
the soft approximation, unlike the collinear approxima- 
tion, is not independent of the line to which the /i end of 
the gluon attaches. Nevertheless, we can apply the same 
soft approximation to all lines in the coUinear-to-light- 
meson jet, since they are all proportional to the same 
light-like vector (up to corrections of order A.Qcu/'mb)- 
For the constituents of each charmonium, we can also 
apply the same soft approximation up to corrections of 
order rric/Q. The reason for this is that, in going from 
the rest frame of the charmonium to the e+e^ CM frame 
or i?-meson rest frame, the momenta of constituents of 
the charmonium undergo boosts that render all of the 
momenta nearly parallel. When we explicitly consider 
those boosts, it follows that, for the charmonium with 
momentum Pi {P2), the plus (minus) components arc 
boosted by a factor of order Q/rric and the minus (plus) 
components are boosted by a factor of order rric/Q, while 
the transverse components are unchanged. Thus, all of 
the momenta of the constituents of the charmonium with 
momentum Pi {P2) are dominated by the plus (minus) 
component, up to corrections of order rric/Q- This result 
holds provided that all of the components are of approxi- 
mately the same size in the charmonium rest frame. This 
is the case for P^, g^, but it is also true for other momenta 
that characterize the charmonium in its rest frame, such 
as a typical potcntial-gluon momentum (which has spa- 
tial components of order mcV and a temporal component 
of order rricv'^) or a typical rest-frame soft-gluon momen- 
tum (which has all components of order rricv). Therefore, 
we use a modified soft approximation for the constituents 
of the charmonium in the jet Ji in which we replace 
5^1^ in the gluon-propagator numerator with k^riii^/k ■ Ui, 
where rii (712) is a vector with unit component in the 
plus (minus) direction and all other components equal 
to zero. This modified soft approximation differs from 
the standard soft approximation for each constituent of 
a quarkonium by terms that are suppressed as rric/Q. 
Hence, it accounts for all soft (logarithmic) singularities, 
up to terms that are suppressed as rric/Q. We can use 
the modified soft approximation, plus the diagrammatic 
Ward identities to decouple the gluons with soft momenta 
from the charmonium jets, relying on the fact that the 
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c and c in each charmonium are in a color-singlet state 
IJ, |l7|, |l8[ . Once the gluons with momenta collinear to 



the jet have been decoupled from H and the gluons with 
soft momenta have been decoupled from the charmonium 
jets, there remain gluons with momenta collinear to the 
jet that begin and end within the jet. In the rest frame 
of the charmonium, these gluons correspond to soft and 
potential gluons. The contributions from regions of mo- 
mentum of order nicV or less that arise from these glu- 
ons (including infrared divergences) can be absorbed into 
nonperturbative NRQCD matrix elements [l^. The as- 
sociated NRQCD short-distance coefficients then contain 
only contributions involving momenta of order rric. The 
NRQCD short-distance coefficients couple to the produc- 
tion process only through the active c and c lines. Thus, 
we have arrived at the factorized form. 

As we have mentioned, the modified soft approxima- 
tion allows us to decouple gluons with soft momentum 
from the charmonium jets, up to corrections that are 
suppressed as nic/Q- Thus, in general, one would ex- 
pect corrections to the factorized form to appear at or- 
der nic/Q. However, in the case of double-charmonium 
production in e^e~ annihilation, S decouples from both 
the coUinear-to-plus charmonium jet and the coUinear- 
to-minus charmonium jet. Each decoupling holds up to 
corrections of order rric/Q, and, so, the overall decou- 
pling holds up to corrections of order m^/s. In perturba- 
tion theory, the factorization-violating corrections may 
be enhanced by logarithms of s/m'^. Furthermore, they 
are infrared divergent. In reality, these infrared diver- 
gences are cut off by nonperturbative effects associated 
with confinement. Our analysis does not determine the 
size of these factorization- violating corrections: it shows 
only that they are proportional to one or two powers of 
rric/Q. The constant of proportionality might be deter- 
mined through experiment and/or lattice calculations. 

At the lowest order in as and v, the decoupling of soft 
gluons from each charmonium is exact. At the lowest 
order in as, the charmonium has only the c and c as 
constituents. At the lowest order in v, one sets qi = 0, 
and the c and c momenta become equal. Then, one can 
apply the same soft approximation for both the c and 
c in the charmonium. Since, at the lowest order in as, 
there is, at most, one soft gluon, the soft decoupling in 
double-charmonium production is exact if either char- 
monium is treated at the lowest order in v. An explicit 
calculation of the one- loop corrections to S'-wave charmo- 
nium production in _B-meson decays at the lowest order 
in V j^] and an explicit calculation of the one-loop cor- 
rections to a[e~^e~ —^ J/ip + 'Xcj] , in which the J/ip is 



20|, confirm our expecta- 



treated at lowest order in 
tion that these corrections should be free of infrared di- 
vergences. P-wave charmonium production in B-meson 



ified in Ref. [J|, we obtain agreement with the calcula- 
tions in Ref. [2l|, which yield only infrared divergences 
that are suppressed as m^/s. As we have discussed, we 
find more generally that the violations of factorization 
are suppressed as rad \fs. At higher orders in v, we ex- 
pect the soft cancellation for each charmonium to hold 
only up to corrections of order md \fs, even in the case 
of S-wave charmonium production. Similarly, in relative 
order and beyond, we expect the soft cancellation for 
each charmonium to hold only up to corrections of order 
rarl \fs, even at the lowest order in v. 
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decays was considered at the one- loop level in Ref. [21 1 
within the factorization framework of Ref. [5| ■ When we 
use the light-cone distribution amplitudes that are spec- 
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